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Abstract 

We write a first order symmetric hyperbolic system coupling the 
Riemann tensor with the derivative of an electromagnetic 2 - form and 
the dynamical acceleration of a perfect relativistic fluid. We determine 
the associated, coupled, Bel - Robinson type energy, and the integral 
equality that it satisfies. 

Dedicated to Guy Boillat with deep esteem and affection. 



1 Introduction. 

The interest of formulating fundamental physical laws as first order symmet- 
ric hyperbolic systems has been stressed, after the foundational work of K. 
O. Friedrichs, by various authors. Particularly important developments in 
the case of systems of conservation laws have been obtained by Lax, Boil- 
lat, Ruggeri and Strumia, Anile, Dafermos. They have used the existence, in 
many physical cases, of an additional convex conserved density, associated to 
an entropy functional. They have obtained symmetric first order hyperbolic 
systems in many cases of physical interest, by introducing auxiliary variables 
similar to Lagrange multipliers. 

In the case of Einstein's gravitation theory, the problem seems more del- 
icate. There is no local density of gravitational energy, nor entropy. The 
effective strength of the gravitational field lies in the Riemann tensor of the 
spacetime metric. Its evolution is governed by the so - called higher order 
equations (Bel 1958, Lichnerowicz 1964), deduced from the Bianchi identi- 
ties. The system satisfied by the trace free part of the Riemann tensor, the 
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Weyl tensor, was some time ago recognized as a linear, first order symmetriz- 
able hyperbolic system (FOSH), with constraints, homogeneous in vacuum. 
See H. Friedrich (1996) and references therein. The evolution equations for 
the Riemann tensor itself, now renamed Bianchi equations, have also been 
written (CB-Yo 1997, 2001) as a FOSH system, made explicit in terms of 
four two- tensors, introduced by Bel (1958), the electric and magnetic grav- 
itational fields and corresponding inductions relative to a Cauchy adapted 
frame. A FOSH system has also been written for the Bianchi equations 
(no longer homogeneous) and the dynamical acceleration of a perfect fluid 
(CB-Yo2002). In this article, dedicated to one of the best specialist of these 
FOSH systems, we couple the Bianchi equations with the equations satisfied 
by the dynamical acceleration of a charged fluid and the derivatives of the 
associated Maxwell field. 

2 Einstein equations with sources. 
2.1 Definitions. 

The Einstein equations, on a 4 - dimensional manifold V, link the Ricci 
tensor of a pseudo Riemannian metric g of Lorentzian signature with the 
stress energy tensor T of sources. These equationsread 

mccz(g)=p = T-^g(tr g T), (2.1) 

that is, in local coordinates, 

R a j3 — Pa/3, with p a/3 = T a/3 — -g a/3 T x . 

The stress energy tensor of an electrically charged, relativistic, perfect 
fluid is the sum of the fluid stress energy tensor: 

t a p = {fi + p)u a u p + pg a(3 

and the Maxwell tensor: 

T cP = F a x F^ - ^g a/3 F x,l F x ^. (2.2) 
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In these formulas u a is the unit kinematical velocity satisfying u a u a — — 1; 
fi,p,S are the specific energy, pressure and entropy. These thermodynamic 
quantities are assumed to be positive and linked by an equation of state: 

A* = MP,S). (2.3) 

The entropy S satisfies a conservation equation 1 along the flow lines: 

u a d a S = 0. (2.4) 

The tensor F is the antisymmetric electromagnetic closed 2 - form. It satisfies 
the Maxwell equations: 

dF = 0, and 5F = J, (2.5) 

where d and 5 are respectively the differential and co differential operator (in 
the metric g), while J is the electric current. We suppose J to be a convection 
current, that is we have: 

J = qu, i.e. J a = qu a , (2.6) 
with q the electric charge density. 

2.2 Fluid equations. 

The conservation laws, consequence of the Einstein equations, read: 

V Q T Q/3 = X7 a (t al3 + r al3 ) = 0. (2.7) 
The Maxwell equations imply that, with J X F^ X called the Lorentz force, 

V a r a/3 = J X F\. (2.8) 
The index / of the fluid, is given by the identity 

/(p,S):=exp / P % , (2.9) 



1 Such an equation can be deduced as usual from a particle number density conservation 
and the Gibbs relation. 
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We introduce the dynamical velocity, tangent vector to the flow lines, 
defined by 

C a := fu a , hence C a C a = -f. (2.10) 

The dynamical velocity incorporates information on the kinematic velocity 
u a and the thermodynamic quantities. The definitions of C and / imply 
that: 

-*VM s f°L = f>g- + f V a j S . (2,1) 
The entropy equation 2.4 and the tensor t a @ become: 

C a d a S = 0, (2.12) 

= (p + p)f- 2 C a C p + g a pp (2.13) 
Therefore the equations 2.7 can be written: 

C^ a [(fi + P )f- 2 C a ] + (fi + p)r 2 C a V a CP + d?p = -J X F? X . (2.14) 

We contract with Cp these equations and we use the equation 2.11 and 
the entropy conservation 2.12 to obtain the following continuity equation 
(under our hypothesis on J, the Lorentz force J A F /3 A is orthogonal to C) 

V a [(fi + p)f- 2 C a }=0, (2.15) 

which can be written: 

V a C* + (/£ - iJ—VaCJ, = 0, (2.16) 

where p! p := |^ is a given function of p and S. 

One deduces from 2.14 and 2.15 the following equations of motion 

C a V a C p + f — - = J —J X F\ (2.17) 

p+p p+p 

Using 2.11 we see that the equations of motion may be written: 

C a {V a Cp - VpC a + (/i + p)-\-C x C x ) 1 *qF Pa } + ~^dpS = 0. (2.18) 
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Remark 2.1 These equations are not independent, their left hand side is 
orthogonal to C. 

In these equations the fluid unknowns are the four components of the 
vector C a , and the scalar S. The specific pressure p is a known function of / 
(i.e. of C a , by 2.10) and of S 

p = p(C a C a ,S), (2.19) 

determined by inverting the relation 2.3. 

We suppose that S is initially constant, it remains then constant during 
the evolution and the term d a S disappears from the equations. It can be 
proved, by adapting the Rendall symmetrization 2 of the perfect uncharged 
fluid, that the equations 2.16 and 2.17 can be written as a first order sym- 
metric hyperbolic system for the dynamical velocity C. The electromagnetic 
field F appears as underivated in these equations. 

3 Bianchi equations. 

In contradistinction with the fluid and electromagnetic sources, the writing 
of the Einstein equations 2.1 as a first order system for the derivatives of 
the metric requires a choice of gauge. It was done straightforwardly by 
Fisher and Marsden using harmonic (now called wave) coordinates. In these 
coordinates the Ricci tensor takes the form of a quasidiagonal, quasilinear 
system of wave equations for the metric g. Though the non linear stability of 
Minkowski spacetime has now been proved (Lindblad and Rodnianski 2005) 
using these coordinates, the fundamental geometrical object signaling the 
existence of a gravitational field remains the Riemann tensor. 

The Bianchi equations satisfied by the Riemann tensor of a metric g, 
solution of the Einstein equations 2.1, are: 

Vo-R/ii.A^ + Vji? fe,A M — V^-Roj.A^ = (3.1) 

and 

— Vo-R°j,A^ — VhR h ;:i,\n = ~^.0M = ~ VoP M A + V^PoA (3.2) 

2 General techniques, using an additional convex entropy function, of Lax 1957, Boillat 
1994, Ruggeri and Strumia 1981 (reported in the book by Anile 1982), could also be used. 
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We recall that a sliced Lorentzian manifold (V, g) is a product M x R 
endowed with a Lorentzian metric g which induces on each M t := M x {t} 
a Riemannian metric gt- The metric g takes in a Cauchy adapted moving 
frame, i.e a frame with its timelike axis orthogonal to the space slices, the 
usual 3+1 form 

ds 2 = -N 2 dt 2 + g tJ {9 1 + (3 l dt){9 3 + fidt). (3.3) 

The gij are the components of g t in a coframe 6 l on M. Particular cases are 
the natural coframe 9 l = dx l , and an orthonormal coframe 9 l = a l -dx^ where 
9ij = ^ij- We denote by g hk the contravariant components of g t , they are, 
in our spacetime frames, also the hk contravariant components of g. The 
derivatives d a are the Pfaff derivatives in the coframe 9° = dt, 9 l = 9 t + (3 l dt, 
that is 

9o = |-/m, a i = 3 i = A>^, (3.4) 
with A the inverse matrix of a whose elements are a*. 

Theorem 3.1 The equations 3.1 and 3.2. are a first order symmetric hyper- 
bolic system for the components of the Riemann tensor in a Cauchy adapted 
frame, with 9 l an orthonormal frame. 

Proof. The principal operator of these equations is diagonal by blocks, 
a block corresponding to the derivatives of the components Ru,\n and Roh,\n. 
for a given pair A, fi, A < /i. Such a block is the following symmetric matrix, 
if 9 l is an orthonormal coframe : 
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d! 
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I 



Each block is symmetric, hence the full principal operator is also sym- 
metric. 

The matrix M l of the coefficients of derivatives is identical to the ma- 
trix M° of the coefficients of the derivatives <9 . It is diagonal, with coefficients 
either 1 or N~ 2 , it is therefore positive definite. ■ 

The energy on a t =constant submanifold M t , constant t, associated to 
the Bianchi equations is the so - called Bel Robinson energy, integral on M t 
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of the positive definite quadratic form M* in the Riemann tensor. If one 
introduces the electric and magnetic fields and inductions space 2- tensors 
associated with the Riemann tensor 3 given by 

rp — pO n = —n n jyhk,lm 

J-'ij — 1 1 i,0j i J-^ij — ^iihk'\ilm rt i 

Hij = 2^ 1 VihkR hk ,oj , Bji = -N 1 7]i hk R j, hk , 

the Bel - Robinson energy density on M t is given by: 

£ g = l -{\E\ 2 + W + \D\ 2 + \B\ 2 ). 

In the case that we are considering of equations with sources, the Bianchi 
equations contain, in addition to the Riemann tensor and its gradient, the 
dynamical acceleration VC and the gradient VF of F, but no derivative of 
these quantities. 



4 Equations for VF. 

The Maxwell equations imply 

5dF + d5F = dJ (4.1) 

that is, using the Ricci identity and some manipulation of indices, the fol- 
lowing semilinear wave equation for F : 

V a V a F 01 = fa, (4.2) 

with 

fa '■= R/3 X F^\ — Ry X Fxp — 2R/3 X n a F Xa + J 7 — V 7 Jp. (4.3) 

It is easy to deduce a symmetric hyperbolic first order system for VF from 
the wave equation satisfied by F. We set 

Fj, a p := V 7 F a/3 , (4.4) 

In a Cauchy adapted frame with orthonormal 9\ the equations to satisfy 
read: 

Ar 2 V F , Q/3 - S^ViF^p = fa (4.5) 

3 rjijk the volume form on M t .. 
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and 

VoF ij0t p — ViF 0:af3 = Roi^Fxp + Roi,i3 X F a \. (4.6) 

The left hand sides of these equations are a symmetric hyperbolic first order 
operator for the four unknowns (Fo >a p, F^ a p) for each pair (a, (3), a < (3. The 
right hand sides contain VC and Riemann(g), but no derivative of these 
quantities. 

The energy density associated to the above system is a positive quadratic 
form, Esjf in VF. It is sometimes called a "superenergy" . 

5 Equations for VC. 

The dynamical acceleration VC satisfies the following equations obtained by 
covariant differentiation of 2.16. and 2.17, and use of the Ricci identities: 

M 7/3 = C a (V Q C 7/3 - W„C ia ) + a 7/3 = (5.1) 

and 

nan '13 

g a ^ a C lP + - l)^-V a C 7/3 + 6 7 = (5.2) 
where we have set 

C 7/3 = V 7 C /3 , (5.3) 



o 7/3 ee C/(C a/3 - C Pa ) + C a C x R-j x + V 7 {(/i + p)- 1 (-C A C A )^C a F /3 J, 

(5.4) 

bj = -R ixC * + V 7 {(^ - 1)— }C Q/3 (5.5) 

The last term in 6 7 is a quadratic form in C Q/ g whose coefficients are functions 
of the C a and S. These functions can be computed by using the identity (we 
use S = constant) 

V 7 /ip ee /4' 2 <9 7 p. 
By the definition of / and the identity 2.4 it holds that 



7a- 
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The equations 5.1. are not independent, because they satisfy the identities 

C p M lP = 0, (5.6) 

the equations 5.1 and 5.2. are not a well posed system. Instead of the 4x4 
equations 4.1 we consider 4 the 4x3 ones: 

M 7 . = Mli - ^M 70 = (5.7) 

The terms in derivatives of C 7 a in these equations can be written in the 
following form: 

C a d a (C^ - ^C 70 ) - (d t - ^d )(C a C 7a ) (5.8) 

Lemma 5.1 The system 5.2,5.7. is equivalent to a FOS (First Order Sym- 
metric) system for C ja with right hand side function of the Riemann tensor 
and VF, not of their derivatives. 

Proof. The system is quasi diagonal by blocks, each block corresponding 
to a given value of the index 7. We will write the principal operator of a block 
by omitting this index. We set 

C 

Ui = Cyi — — C 7 o, Uq = C a C ia (5-9) 

and we define the differential operators d a as follows: 

d = C a d a , d t = d t -^-d (5.10) 

The principal terms (derivatives of C 7Q ) in the equations 5.7 with index 7 
are 

doUi-diUo. (5.11) 

We have by inverting 5.9: 

_ _ C (U - C l Uj) 
^70= 



4 An analogous procedure is used for the symmetrization of the Eulcr equations in K.O. 
Friedrichs 1969 and in Rcndall 1992. 
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a 



Q(U - CWj) 



C X C X 

The principal terms of 5.2. read, using the above formulae 



v' p c a d a u 



)d i U j 



c o C i 



d U, t 



C X C\ vy C X C X J J C X C X 
We introduce the positive definite (if C is timelike) quadratic form 



(5.12) 



Then we find that 



i«J — M 



c x c x 

c o C i 



(5.13) 



J C C X C\ 
The principal terms 5.12 are therefore 



~c x cT 



(5.14) 



The matrix of the coefficients of the derivatives d a in the equations deduced 
from the system 5.2, 5.7 is 



/- 


MP p 


-o 1 

do 


-d 2 




-d 3 \ 











do 





V 


-h 








Bo ) 



This matrix is symmetrized by taking the product with the 4x4 matrix 

1 

9ij 
Hyperbolicity. 

In the case we are considering the matrix M° is diagonal, with positive 
elements if /i' p > and C a is timelike. The principal matrix written with 
natural coordinates is also symmetrizable. The corresponding matrix M* is 
not diagonal, and it is not obvious that it is positive definite. In fact it will 
be so only if fj,' p > 1. We will prove the following lemma. 



Lemma 5.2 The system 5.2, 5.7 is FOSH if n' > 1 and C is timelike. 
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Proof. It is simpler to compute directly the energy density for the 
considered system: its positivity is equivalent to the positivity of the matrix 
M l . Multiplying 5.2. and 5.7 respectively by Uq and g^Uj gives equations of 
the form 



2 p 



-g t3 U d i U j = U Q (5.15) 



gPUjdaUi - g t3 UjdiU Q = g^U^ (5.16) 

where the $ a contain only non differentiated terms in U, Riemann and VF. 
We add these two equations, replace the operators d by the operators d and 
carry out some manipulations using the expression for g lJ and the Leibniz 
rule. We obtain that 

d £ m + V i ?V m = Q m (5.17) 

The function Q m is a quadratic form in C ia , F a ^ y and the Riemann tensor, 
while S m is the energy density on M t of the dynamical acceleration VC. It 
is the quadratic form given by: 

£ m = f- 2 W P - l)^o + (U - C*U t ) 2 ] + g'UliUj (5.18) 
It is positive definite if fi' > 1 and C a is timelike. ■ 

Remark 5.3 The system is hyperbolic in the sense of Leray if n' p > 0, but 
the submanifolds x° = constant are 'spacelike' with respect to the fluid wave 
cone only if the fluid sound speed is less than the speed of light, i.e., /a' > 1. 

The quantity £ m is called the fluid "acceleration energy" density. 



6 Coupled system. 

The previous results give, for Rieman(g), VF and VC, when g, a lorentzian 
metric, F, a 2 form and C, a timelike vector, are known, a first order sym- 
metric system, quasi diagonal by blocks, hyperbolic if /i' > 1. 

By choosing a densitized lapse (time wave gauge), as in [CB-Yo97] and 
[CB-YoOl] one can obtain a full symmetric first order system containing, in 
addition, g and its connection, F and C. 
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The superenergy density 8 for the full system is the sum of the Bel - 
Robinson energy density of the gravitational field, the energy density of VF, 
supernergy of F, and the energy density of the dynamical acceleration: 



Using the expression of do and the mean extrinsic curvature r = g^Kij 
of the space slices S t , whose volume element we denote by fig t , we obtain 
by integration an integral equality whose right hand side couples all these 
superenergies: 



where Q is a quadratic form in Riemann(g), VF and VC, which could be 
estimated in terms of £ if the other unknowns, C, g and its connection, were 
estimated. 
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